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ANALYSIS  OF  ACOUSTIC  IMPEDANCE  DATA 


Abstract:  Procedures  and  applicable  charts  for  the  analysis  of 
acoustic  impedance  data  obtained  from  the  USRL  pulse  tube  are 
presented.  Formulas  and  tables  for  rapid  determination  of  the 
bilinear  transformation  (l  +  r)/(l  -  r)  are  shown. 


INTRODUCTION 

The  USRL  has  recently  completed  development  of  a  system  for  measuring 
the  complex  acoustic  impedance  and  the  complex  propagation  constant  of 
acoustical  materials  for  underwater  sound.  Basically,  the  system  is  an 
acoustic  transmission  line  180  cm  long  with  a  reversible  transducer  at 
one  end  and  a  sample  of  the  acoustical  material  of  interest  at  the  oppo¬ 
site  end.  The  operating  frequency  range  is  2.5  to  approximately  10  kHz. 
Hydrostatic  pressure  is  variable  to  10,000  psi. 

For  a  sample  inserted  in  the  180-cm-long  acoustic  transmission  lino 
of  the  system,  two  quantities  are  measured  at  each  operating  frequency 
and  hydrostatic  pressure.  These  quantities  are  the  magnitude-  |r|  and  the 
phase  angle  0  of  the  complex  reflection  coefficient  r  at  the  front  bound¬ 
ary  of  the  sample. 

This  report  describes  the  details  of  the  analysis  to  be  performed  on 
the  measured  data  |r|  and  e  to  calculate  the  two  constants  that  describe 
the  acoustical  behavior  of  the  sample  at  each  discrete  frequency  and 
hydrostatic  pressure.  These  are  the  attenuation  constant  a  (usually 
expressed  in  dB/m)  and  the  speed  of  sound  propagation  c  in  the  sample. 

The  basic  equation  that  describes  the  wave  propagation  in  a  sample 
terminated  by  a  rigid  boundary  is  (see  Appendix  for  derivation) 

(Zin/ZT)  =  (Zq/Z-j.)  coth  vd,  (1) 

where  Z. is  the  normalized  acoustic  impedance  at  the  water-to-sample 
boundary,  Z^Z^  is  the  normalized  characteristic  impedance  of  the  sample, 
Zj  =  pTcT  is  the  characteristic  impedance  of  water,  pT  is  the  density  of 
water,  c^  is  the  speed  of  sound  in  water,  y  =  a  +  jk  is  the  complex  pro¬ 
pagation  constant  of  the  sample,  k  =  2Trf/c  is  the  wavenumber  in  the  sam¬ 
ple,  d  is  the  length  of  the  sample,  and  f  is  the  frequency. 


The  ratio  ZVZ_  in  Eq.  (l)  can  be  represented  by  (see  Appendix  for 
derivation) 


ZT  pT  a  +  jk 


where  k^  =  2nf/c^.  is  the  wavenumber  in  water  and  p  is  the  density  of  the 

sample.  The  measured  complex  reflection  coefficient  r  for  the  sound 
pressure  represents  an  indirect  measurement  of  the  water- to-sampls  bound¬ 
ary  impedance  according  to  the  relation 

(Zin/ZT)  =  (l  +  r)/(l  -  r).  (3) 

For  brevity,  the  symbol  Z  will  be  used  henceforth  to  represent  Z^/Zy. 

When  the  right-hand  terms  of  Eqs.  (2)  and  (3)  are  substituted  into 
Eq.  (l),  the  following  relation  is  obtained* 


1  +  r  p  k-d 

- =  j  —  -  coth  (a  +  jk)d.  (4) 

1  -  r  pT  (a  +  jk)d 

In  addition  to  r,  the  following  quantities  in  this  equation  are  measured: 

p/p  ,  the  ratio  of  the  density  of  the  sample  to  that  of  water;  and  k  d, 

T  •  ^ 

the  acoustic  length  in  radians  of  a  column  of  water  whose  length  is  equal 

to  that  of  the  sample.  Thus,  when  all  of  the  measured  quantities  are 

written  on  the  left-hand  side,  Eq.  (4)  becomes 


1  1  +  r  coth  (o'  +  jk)d 

- - =  j  - 

(p/p  )k  d  1  -  r  (a  +  jk)d 


Because  Eq.  (5)  is  a  transcendental  equation,  it  cannot  be  solved  for  the 
unknown  values  of  o’  and  k  by  algebraic  methods.  The  procedure  is  to  use 
a  chart  or  family  of  charts  presenting  the  value  of  the  complex  function 
F  for  assumed  values  of  the  complex  independent  variable — in  this  case, 
ad  +  jkd.  Thus  F  is  defined  by 


F  = 


coth  (a  +  jk)d 
(a  +  jk)d 


(6a) 


and  is  related  to  Z  by 


(p/pT)kTd 


(6b) 


and 


Imag  F  Imag  Z 

tan  <f>  = - = - 

Real  F  Real  Z 


(6c) 


Computations  were  made  on  an  IBM  1620  computer  to  evaluate  the  func¬ 
tion  F  by  Eq.  (6a)  for  many  values  of  2kd  and  Q,  where  Q  =  kd/2cyd.  In 
these  computations  2kd  is  in  radians  and  Q  is  derived  from  a  value  of  2od 
in  nepers.  The  computed  values  |F|  and  (imag  F)/(Real  F),  when  plotted 
logarithmically,  yield  the  charts  presented  in  Figs.  l(a)  through  1(e). 

The  user  enters  the  appropriate  chart  through  the  rectangular  coordinates 
with  the  measured  value,  the  left-hand  side  of  Eq.  (5).  The  corresponding 
values  of  2kd  and  Q  are  read  from  the  curvilinear  coordinates.  From  the 
values  of  2kd  and  Q,  the  constants  are  computed  by 

2kd  8.686 

a  = - -  - dB/m,  (7a) 

Q  4d 

c  =  4rcdf/2kd  m/sec.  (7b) 


RAPID  METHOD  FOR  THE  TRANSFORMATION  (l  +  r)/(l  -  r) 

The  left-hand  side  of  Eq.  (5)  shows  that  the  measured  data  (reflec¬ 
tion  coefficient  r)  have  to  be  transformed  by  the  expression  (1  +  r )/( 1  ~r). 
The  Smith  chart  of  transmission  line  theory  is  frequently  used  for  this 
transformation.  In  our  case,  we  want  the  impedance  (1  +  r)/(l  -  r)  in 

the  polar  form  Z  =  |Z|e^.  The  appropriate  chart  to  obtain  impedance  in 
polar  form  is  the  modified  Smith  chart  described  by  Beranek  [1]. 

The  result  of  this  r-to-Z  transformation  is  to  be  used,  after  division 
by  (p/p^)k^d  as  in  Eq.  (6b),  in  the  charts  of  Fig.  1.  The  successive  use 

of  two  separate  charts  (the  modified  Smith  chart  and  the  charts  of  Fig.  l) 
is  inadvisable  from  the  point  of  view  of  final  accuracy.  For  this  reason 
the  r-to-Z  transformation  should  be  performed  numerically.  From  the 
relation 

Z  =  (1  +  r)/(l  -  r),  (8) 


the  expressions  for  jz|  and  are 

fl  +  2  |r  |  cos  6  -»•  |  r  |  ^1^ 


Z|  = 


2lrlcos  0  +  I r I 


2  ' 


and 


2 Ir | sin  9 

tan  <f>  =  - .  (9b) 

1  -  |r|2 

Computations  by  Eqs.  (9a)  and  (9b)  are  lengthy  and  laborious.  It  is 
possible,  however,  to  modify  the  equations  for  simpler  computations.  The 
introduction  of  the  auxiliary  angle  R,  defined  by  tan  -gR  =  |r|,  simplifies 
Eqs.  (9a)  and  (9b)  into 


.0 


IMAGINARY  F 
REAL  F 


I  I  I  l  I  I  I 
-0.1 


Fig.  1(e).  Chart  of  the  computed  function 


coth  ( oi  +  jk)d 

p  =  j  - • 

( a  +  jk)d 


1 

7 


> 


(10a) 


1 

+  sin 

R 

cos 

0 

1 

-  sin 

R 

cos 

0. 

and 


tan  0  =  tan  R  sin  8. 


(10b) 


The  auxiliary  angle  R  is  always  in  the  first  quadrant  (because  |r|  £  l) 
and  accordingly  both  sin  R  and  tan  R  in  Eqs.  (10a)  and  (10b)  are  inher¬ 
ently  positive  in  all  cases.  Thus,  by  Eq.  (10b),  tan  0  obtains  the  sign 
of  sin  0;  this  relation  together  with  the  condition  that  0  is  restricted 

to  the  first  or  to  the  fourth  quadrant  (because  |Z|ei0  physically  is  a 
real  impedance)  defines  the  specific  quadrant  of  0  unambiguously.  When  9 
is  in  the  first  or  in  the  second  quadrant,  0  is  in  the  first  quadrant; 
when  6  is  in  the  third  or  in  the  fourth  quadrant,  0  is  in  the  fourth 
quadrant. 


By  regarding  |sin  R  cos  0|  as  the  cosine  of  another  auxiliary  angle  S, 
(S  restricted  to  the  first  quadrant  by  definition),  Eq.  (10a)  becomes 


=  cot  -g-S,  if  sin  R  cos  0  is  positive,  ( 1 1  a ) 


[j  +  cos  S 


=  tan  ^S,  if  sin  R  cos  0  is  negative.  (lib) 


It  follows  from  Eqs.  ( 1 1  a )  and  (lib)  that  |Z|  >  1  when  9  is  in  the  first, 
or  in  the  fourth  quadrants,  and  |Z|  <  1  when  0  is  in  the  second  or  in  the 
third  quadrants. 


Tables  I  and  II  are  intended  for  use  with  Eqs.  ( 1 0)  and  (ll).  Since 
the  rectangular  coordinates  of  the  Fig.  1  charts  are  logarithmic,  the 
tables  are  also  logarithmic  to  facilitate  the  plotting  of  measured  data 
points.  In  our  measurements  |r|  was  measured  to  the  nearest  0.1  dB;  for 
this  reason  Table  I  shows  values  of  I r |  in  the  form  20  log  |r|  in  incre¬ 
ments  of  0.1  dB  from  -0.1  to  -6.0  dB. 


A  representative  example  extracted  from  a  series  of  measurements  with 
the  recently  developed  acoustic-impedance-measuring  system  will  show 
details  of  the  r-to-Z  computations  by  the  method  described. 

The  measured  quantities  are  20  log  |r[  =  -2.6  dB  and  0  =  105°-  Then, 


from  trig  tables: 

log  sin  0  =  9«9849(pos) 

log  cos  0  =  9*4130(neg) 

From  Table  I: 

loq  tan  R  =  0.5174(oos) 

loq  sin  R  =  9.9808(oos) 

Add: 

log  tan  0  =  O.5023(pos) 

log  cos  S  =  9.3938(neg) 

From  Table  II: 

log  |Z|  =  9.8901 

Result: 

|Z|  =  0.776,  tan  0  =  3-18. 
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For  clarification  it  is  pointed  out  that  (pos)  and  (neg)  are  used  above 
to  keep  track  of  the  sign  of  the  quantities.  For  example,  log  sin  6  = 
9.9849(pos)  indicates  that  sin  0  is  positive  and  log  cos  8  =  9.4130  (neg) 
is  used  to  indicate  that  cos  9  is  negative. 

This  completes  the  r-to-Z  computation.  To  enter  the  chart,  the  value 
of  (p/py)kyd,  where  kyd  is  in  radians,  is  needed  to  compute  |Fj  as  in 

Eq.  (6b).  In  the  example  being  discussed,  p/py  =  1.41  and 

k^,d  =  2.48  radians,  so  |F|  =  0.222.  The  chart  (Fig.  l(b)  is  used  in  this 

example)  is  entered  through  the  rectangular  coordinates  with  the  values 
tan  0  =  3*18  and  |F  |  =  0.222;  Eq.  (6c)  points  out  that  tan  0  is 
(imag  F)/(Real  F). 

For  the  precise  location  of  the  point  that  represents  tan  0  =  3*18 
and  | F J  =  0.222  it  must  be  kept  in  mind  that  the  charts  are  logarithmic, 
hence  logarithmic  interpolation  must  be  used  between  the  major  scale 
divisions  of  the  charts. 

It  is  equally  valid  to  consider  the  charts  as  linear  plots  of  log  |F| 
and  log  tan  0.  For  this  viewpoint  the  reader  needs  to  visualize  that 
values  of  log  |F|  and  log  tan  0  can  be  plotted  directly  on  the  chart  in  a 
linear  fashion.  Still  for  the  same  example,  log  tan  0  =  O.5023(pos)  and 
log  j F |  =  9.3465(pos);  the  chart  of  Fig.  l(b)  is  entered  on  the  abscissa 
between  1  and  10  (corresponding  to  the  characteristic  0  of  log  tan  0)  at 
0.502  of  the  distance  between  1  and  10.  For  the  ordinate  on  the  chart, 
the  entering  point  is  between  0.1  and  1  (corresponding  to  the  characteris¬ 
tic  9  of  log  [F [ )  at  0.346  of  the  distance  between  0.1  and  1.  On  the 
curvilinear  coordinates,  the  plotted  point  locates  the  values  2kd  =  3*95 
and  Q  =  7.  To  complete  the  computation,  Eqs.  (7)  are  used  to  find  a  and 
c  (d  =  0.100  m,  f  =  6.0  kHz,  in  this  example), 

(3.95) (8.686) 

a  =  -  =  12  dB/m, 

7(0.400) 

4n(0.100)f 

c  = - =  1910  m/sec. 

3.95 

A  few  final  remarks  are  probably  in  order  concerning  Fig.  1.  The 
charts  could  have  been  plotted  in  any  one  of  several  other  forms.  The 
form  presented  was  found  most  expedient  for  use  with  our  measurements. 

A  possible  source  of  confusion  to  the  reader  may  be  the  presence  of  linear 
scale  divisions  on  an  otherwise  logarithmic  chart.  The  example  above  has 
attempted  to  show  that  the  charts  may  be  regarded  as  logarithmic  for 
tan  0  and  | F J  or  as  linear  for  log  tan  0  and  log  |F|. 

The  inquisitive  reader  is  referred  to  the  works  of  other  authors 
[4-10]  for  discussions  of  the  mapping  of  the  (coth  w)/w  and  (tanh  w)/w 
functions  in  connection  with  a  variety  of  applications  in  airborne  acous¬ 
tics  and  in  microwave  dielectric  measurements. 
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IP1 


R  and 
f  1  r  | 


log 
sin  R 


log 
sin  R 


■0.1  1.9386  0.0000 

-0.2  1.6377  9.9999 
•0.3  1.4616  9.9997 
-0.4  1.3366  9.9995 


1.2396 

1.1603 

1.0932 

1.0351 

0.9838 


9.9993 

9.9990 

9.9986 

9.9982 

9.9977 

9.9971 

9.9965 

9.9959 

9.9952 

9.9944 

9.9936 

9.9927 

9.9917 

9.9907 

9.9897 

9.9886 

9.9874 

9.9862 

9.9850 

9.9836 

9.9822 

9.9808 

9.9794 

9.9778 

9.9762 


0.4531 

9.9746 

0.4383 

9.9729 

0.4239 

9.9712 

0.4099 

9.9694 

0.3963 

9.9675 

0.3831 

9.9657 

0.3701 

9.9637 

0.3576 

9.9617 

0.3453 

9.9597 

0.3333 

9.9576 

0.3215 

9.9555 

0.3100 

9.9533 

0.2988 

9.9511 

0.2877 

9.9488 

0.2769 

9.9465 

0.2663 

9.9442 

0.2559 

9.9418 

0.2457 

9.9393 

0.2357 

9.9368 

0.2258 

9.9343 

0.2161 

9.9317 

0.2066 

9.9291 

0.1972 

9.9264 

0.1879 

9.9237 

0.1788 

9.9210 

0.1698 

9.9182 

0.1609 

9.9154 

0.1522 

9.9125 

0.1436 

9.9096 

0.1350 

9.9067 

C.1266 

9.9037 

mm mm 


Ij55£gyg25=jj 


Table  II 


Tabulation  for  Finding  (Z 


log 
cos  S 


log 

cot  -g-S 


log 
cos  S 


0 
2 
4 

8.6000 

8.8000 

9.0000 

9.1000 

9.2000 

9.3000 

9.3500 

9.4000 

9.4500 

9.5000 

9.5500 

9-6000 

9.6200 


9.7000 

9.7200 

9.7400 

9.7600 

9.7800 

9.8000 

9.8100 

9.8200 

9.8300 

9.8400 

9.8500 

9.8600 

9.8700 

9.8800 

9.8900 


0.0044 

9.9956 

0.0069 

9.9931 

0.0109 

9.9891 

0.0174 

9.9826 

0.0273 

9.9727 

0.0436 

9.9564 

0.0550 

9.9450 

0.0694 

9.9306 

0.0878 

9.9122 

0.0989 

9.9011 

0.1115 

9.8885 

0.1258 

9.8742 

0.1422 

9.8578 

0.1611 

9.8389 

0.1830 

9.8170 

0.1928 

9.8072 

0.2032 

9.7968 

0.2144 

9.7856 

0.2264 

9.7736 

0.2392 

9.7608 

0.2532 

9.7468 

0.2683 

9.7317 

0.2847 

9.7153 

0.3028 

9.6972 

0.3227 

9.6773 

0.3334 

9.6666 

0.3448 

9.6552 

0.3569 

9.6431 

0.3698 

9.6302 

0.3835 

9.6165 

0.3982 

9.6018 

0.4140 

9.5860 

0.4312 

9.5688 

0.4499 

9.5501 

9.9000  0.4704 
9.9050  0.4814 


9.9100 

9.9150 

9.9200 

9.9250 

9.9300 

9.9350 

9.9400 

9.9450 

9.9500 

9.9550 

9.9600 

9.9620 

9.9640 

9.9660 

9.9680 

9.9700 

9.9720 

9.9740 

9.9760 

9.9780 

9.9800 

9.9810 

9.9820 

9.9830 

9.9840 


0.4931 

0.5054 

0.5185 

0.5324 

0.5473 

0.5634 

0.5807 

0.5995 

0.6202 

0.6430 

0.6685 

0.6797 

0.6914 

0.7038 

0.7169 

0.7309 

0.7459 

0.7620 

0.7794 

0.7982. 


0.8674 

0.8814 

0.8963 

0.9124 

0.9299 

0.9488 


9.5296 

9-5186 

9.5069 

9.4946 

9.4815 

9.4676 

9.4527 

9*4366 

9.4193 

9.4005 

9.3798 

9.3570 

9.3315 

9.3203 

9.3086 

9.2962 

9.2831 

9.2691 

9.2541 

9.2380 

9.2206 

9.2018 

9.1812 

9.1699 

9.1582 

9.1458 

9.1326 

9.1186 

9.1037 

9.0876 

9.0701 

9.0512 


log 

log 

log 

cos  S 

cot  -§-S 

tan  -g-S 

9.9900 

0.9694 

9.0306 

9.9905 

0.9806 

9.0104 

9.9910 

0.9922 

9.0078 

9.9915 

1.0047 

8.9953 

9.9920 

1.0179 

8.9821 

9.9925 

1.0319 

8.9681 

9.9930 

1.0469 

8.9531 

9.9935 

1.0630 

8.9370 

9.9940 

1 .0803 

8.9197 

9.9945 

1 .0992 

8.9008 

9.9950 

1.1199 

8.8801 

9.9952 

1.1288 

8.8712 

9.9954 

1.1380 

8.8620 

9.9956 

1.1477 

8.8523 

9.9958 

1.1578 

8.8422 

9.9960 

1.1684 

8.8316 

9.9962 

1.1795 

8.8205 

9.9964 

1.1913 

8.8087 

9.9966 

1.2037 

8.7963 

9.9968 

1 .2168 

8.7832 

9.9970 

1.2308 

8.7692 

9.9972 

1.2458 

8.7542 

9.9974 

1.2619 

8.7381 

9.9976 

1.2793 

8.7207 

9.9978 

1.2982 

8.7018 

9.9980 

1.3189 

8.6811 

9.9982 

1.341 8 

8.6582 

9.9984 

1.3674 

8.6326 

9.9986 

1 .3953 

8.6047 

9.9988 

1 .4298 

8.5702 

9.9990 

1 .4694 

8. 5306 

9.9991 

1 .4923 

8.5077 

9.9992 

1.5178 

8.4822 

9.9993 

1.5469 

8.4531 

9.9994 

1.5803 

8.4197 

9.9995 

1.6199 

8.3801 

9.9996 

1 . 6684 

8.3316 

9.9997 

1.7308 

8.2692 

9.9998 

1.8188 

8.1812 

9.9999 

1 .9694 

8. 0306 

APPENDIX 


From  transmission  line  theory,  the  equation  for  the  input  impedance 
of  a  line  section  terminated  by  an  open-circuit  (analogous  to  a  rigid 

boundary- -that  is,  high  impedance)  is  [2]: 

Z.  =  Z  coth  yd,  (A1 ) 

in  0  T 

where  Z^  and  y  are  the  characteristic  impedance  and  propagation  constant, 

respectively,  of  the  line  section,  and  d  is  its  length.  When  this  line 
section  is  used  as  the  termination  impedance  of  another  transmission  line 
whose  characteristic  impedance  is  Zp  then  Eq.  (l)  describes  the  termina¬ 
tion  impedance  normalized  to  the  characteristic  impedance  Z . 

Equation  (2)  car-  be  derived  from  Eq.  (B4)  in  McSkimin  [3],  which  is 


z0  =  jpu/(<*  +  jk). 

When  both  sides  of  the  equation  are  divided  by  Zp  recalling  that 
Z =  pTCp  Eq.  (A2)  becomes 


0  _ 


ZT  PTcT(a  +  jk) 

and,  since  to/c  is  kp  Eq.  (2)  is  verified. 


(A2) 


(A3) 
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